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We investigate predictions of a minimal realistic non-supersymmetric SU(5) grand unified 
theory. To accomplish unification and generate neutrino mass we introduce one extra Higgs 
■ representation — a 15 of SU(5) — to the particle content of the minimal Georgi-Glashow sce- 

nario. Generic prediction of this setup is a set of rather light scalar leptoquarks. In the case 

>' 

\q | of the most natural implementation of the type II see-saw mechanism their mass is in the 
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phenomenologically interesting region (O(10 2 -10 3 ) GeV). As such, our scenario has a poten- 
tial to be tested at the next generation of collider experiments, particularly at the Large 
■ Hadron Collider (LHC) at CERN. The presence of the 15 generates additional contributions 

O 

gauge d = 6 ones. We exhaustively study both and show that the scenario is not excluded 
by current experimental bounds on nucleon lifetimes. 
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jj ■ I. INTRODUCTION 

The most predictive grand unified theory (GUT) based on an SU(5) gauge symmetry is a 
minimal non-supersymmetric model of Georgi and Glashow 1] (GG). However, the failure to 
accommodate experimentally observed fermion masses and mixing and to unify electroweak and 
strong forces decisively rule it out. Nevertheless, the main features of the underlying theory, e.g, 
partial matter unification and one-step symmetry breaking, are so appealing that there has been a 
number of proposals to enlarge its structure by adding more representations to have the theory in 
agreement with experimental data. Since the number of possible extensions is large it is important 
to answer the following question. What is the minimal number of extra particles that renders 
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an SU(5) gauge theory realistic? Such an extension of the GG model with the smallest possible 
number of extra particles that furnish full SU(5) representation(s) has all prerequisites to be the 
most predictive one. Thus, if there is a definite answer to the first question it is important to 
ask the second one: What are the possible experimental signatures and associated uncertainties of 
such a minimal setup? If uncertainties of the minimal extension are significant the same is even 
more true of a more complicated structure unless additional assumptions are imposed. We address 
both questions in great detail and present truly minimal, i.e., minimal in terms of number of fields, 
realistic non-supersymmetric SU(5) scenario. 

As we demonstrate later, unification, in the minimal scenario, points towards existence of light 
scalar leptoquarks which could generate very rich phenomenological signatures. These are looked 
for in the direct search experiments as well as in the experiments looking for rare processes. Their 
presence generates novel proton decay contributions which can be very important. Moreover, 
since in our scheme their coupling to matter is through the Majorana neutrino Yukawas, their 
observation might even allow measurements of and/or provide constraints on neutrino Yukawa 
coupling entries. This makes our scenario extremely attractive. At the same time, rather low scale 
of vector leptoquarks that is inherent in non-supersymmetric theories exposes our scenario to the 
tests via nucleon lifetime measurements. We investigate all relevant experimental signatures of the 
scenario, including its status with respect to the present bounds on nucleon lifetime. 

In the next section we define our framework. In Section ITTT1 we discuss how it is possible to get 
gauge couplings unification in agreement with low energy data. Then, in Section IIVI we discuss 
possible experimental signatures of the minimal scenario. We conclude in Section Appendix lAl 
contains relevant details and notation of the minimal non-supersymmetric realistic SU(5) we refer 
to throughout the manuscript. The origin of theoretical bounds on the familiar gauge d = 6 proton 
decay operators is critically analyzed in Appendix[Bl Appendix[0 contains details on the two- loop 
running of the gauge couplings that is presented for completeness of our work. 

II. A MINIMAL REALISTIC NON-SUPERSYMMETRIC SU(5) SCENARIO 

In order to motivate the minimal SU(5) grand unifying theory (GUT), where we define such 
a theory to be the one with the smallest possible particle content that renders it realistic, we 
first revisit the GG model jl| and discuss its shortcomings. Only then do we present the minimal 
realistic scenario and investigate its experimental signatures and related uncertainties. 

The GG model fails from the phenomenological point of view for a number of reasons: 
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1. It does not incorporate massive neutrinos; 

2. It yields charged fermion mass ratios in gross violation of experimentally observed values; 

3. It cannot account for the gauge coupling unification. (For one of the first rulings on 
(non)unification see for example Ref. |2j.) 

The first flaw is easy to fix; one either introduces three right-handed neutrinos — singlets of the 

n 

Standard Model (SM) — to use a type I see-saw mechanism 3] to generate their mass or adds a Higgs 
field — a 15 of SU(5) — to generate neutrino mass through the so-called type II see-saw 4,1^. One 
might also use the combination of the two. The fourth option to use the Planck mass suppressed 
higher-dimensional operators jf], Q] does not look promising since it generates too small scale for 
neutrino mass to explain the atmospheric and solar neutrino data. Nevertheless, it might still 
play an important role We focus our attention on the second option, i.e., addition of 15, in 
view of the fact that the right-handed neutrinos, being singlets of the SM, do not contribute to 
the running. Hence, their mass scale cannot be sufficiently well determined or constrained unless 
additional assumptions are introduced. 

The second flaw can be fixed by either introducing the higher-dimensional operators in the 

q , n 

Yukawa sector |9( or resorting to a more complicated Higgs sector a la Georgi and Jarlskog |1(J |. 
The former approach introduces a lot more parameters into the model (see for example |ll|]) but 
unlike in the neutrino case these operators might have just a right strength to modify "bad" mass 
predictions for the charged fermions. In order to keep as minimal as possible the number of particles 
we opt for the scenario with the non-renormalizable terms. 

The third flaw requires presence of additional non-trivial split representations besides those of 
the GG SU(5) model. It can thus be fixed in conjunction with the first and second one. For 
example, introduction of an extra 45 of Higgs to fix mass ratios of charged fermions a la Georgi 
and Jarlskog 10] allows one to achieve unification and, at the same time, raise the scale relevant for 
proton decay. (For the studies on the influence of an extra 45 on the running and other predictions 



see for example 
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13j.) We will see that the addition of one extra 15 of Higgs plays a crucial 
role in achieving the unification in our case. 

So, what we have in mind as the minimal realistic SU(5) model is the GG model supplemented 
by the 15 of Higgs to generate neutrino mass and which incorporates non-renormalizable effects to 
fix the Yukawa sector of charged fermions. We analyze exclusively non-supersymmetric scenario 
for the following three reasons: first, this guarantees the minimality of the number of fields; second, 
there are no problems with the d = 4 and d = 5 proton decay operators; third, since the grand 
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unifying scale is lower than in supersymmetric scenario the setup could possibly be verified or 
excluded in the next generation of proton decay experiments. 

III. UNIFICATION OF GAUGE COUPLINGS 

The main prediction, besides the proton decay, of any GUT is the unification of the strong and 
electroweak forces. We thus show that it is possible to achieve gauge coupling unification in a 
consistent way in our scenario. 

At the one-loop level the running of gauge couplings is given by 

_ii -1.1,, M GUT , . 

a i \m z = a GUT + —hi In (1) 

where i = 1,2,3 for U(l), SU(2), and SU(3), respectively, bi are the appropriate one-loop coeffi- 
cients |3| and acuT = Qgut/^ 71 ) represents the gauge coupling at the unifying scale Mqut- The 
SM coefficients for the case of n light Higgs doublet fields are: 

40 n 20 n , 

h = 1 , b 2 = 1 — , b 3 = -7. (2) 

10 10' 6 6' K J 

Even though the SM coefficients do not generate unification in both the n = 1 and n = 2 case for 

any value of acur and Mqut that is not an issue since the SM does not predict the gauge coupling 

unification in the first place. On the other hand, a GUT, which does predict one, automatically 

introduces a number of additional particles with respect to the SM case that, if light enough, can 

change the outcome of the SM running. This change is easily incorporated if one replaces bi in 

Eq. ffl with the effective one-loop coefficients Bi defined by 



where bn are the one-loop coefficients of any additional particle I of mass Mi {Mz < Mi < Mqjjt)- 
Basically, given a particle content of the GUT and Eqs. (0) and © we can investigate if the 
unification is possible. 

Following Giveon et al jl^ |. Eqs. (pQ) can be further rewritten in a more suitable form in terms 
of differences in the effective coefficients Bij{= Bi — Bj) and low energy observables. They find 
two relations that hold at Mz- 

B 23 5 sin 2 9 W - a em /a s 



B 12 8 3/8 - sin 



2 a 

w 



(4a) 



ln M GUT = 16tt 3/8 - sin 2 6 W 
Mz 5a em B\2 



Adopting the following experimental values at Mz in the MS scheme sin 2 9 W = 0.23120 ± 
0.00015, a'^ = 127.906 ± 0.019 and a s = 0.1187 ± 0.002, we obtain 

= 0.719 ±0.005, (5a) 

-D12 

Mqut 184.9 ±0.2 

= — (5b) 

Last two equations allow us to constrain the mass spectrum of additional particles that leads to 
an exact unification at Mqut- (In what follows we consistently use central values presented in 
Eqs. (JSJ) unless specified otherwise. The inclusion of the two-loop effects and threshold corrections 
is addressed in detail in Appendix 0) 

The fact that the SM with one (two) Higgs doublet (s) cannot yield unification is now more 
transparent in light of Eq. (|5a|). Namely, the resulting SM ratio is simply too small (-B23/-B12 = 0.53 
for n = 1) to satisfy equality in Eq. (|5a|). What is needed is one or more particles that are relatively 
light and with suitable 6« coefficients that can increase the value of the -B23/-B12 ratio. The most 
efficient enhancement is realized by a field that increases -B23 and decreases Byi simultaneously. 
For example, light Higgs doublet is such a field (see the \l/£)(c 5h) coefficients in Table [IJ and it 
takes at least eight of them, at the one-loop level, to bring -B23/-B12 in accord with experiments. 
Other fields that could generate the same type of improvement in our scenario are light 24h), 
<3? a (C 15h) and ^(C 15#). By coefficients of all the particles in our scenario are presented in 
Table [I] and the relevant notation is set in Appendix 

TABLE I: Contributions to the Bi 3 coefficients. The mass of the SM Higgs doublet is taken to be at Mz- 





Higgsless SM 






V 




s 3 




$6 




B23 

B\2 


11 

3 

22 
3 


1 

6 

1 

15 


T5 r *T 


-fry 
-7r v 











IS'*. 



The improvement can also be due to the field that lowers B12 only or lowers B12 at sufficiently 
faster rate than i?23- Looking at Table [J we see that the superheavy gauge fields V comprising X 
and Y gauge bosons and their conjugate partners can accomplish the latter. (Note that the gauge 
contribution improves unification at the one-loop level only if n 7^ and the improvement grows 

, n|. This is because the B23/B12 ratio of the Higgsless SM coefficients 



with the increase of n 
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is the same as for the corresponding ratio of V coefficients.) But, their contribution to running 
has to be subdominant; otherwise one runs into conflict with the experimental data on nucleon 
lifetimes. 
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All in all, the fields capable of improving unification in our minimal SU(5) grand unified scenario 
are ^?d, S3, <£ a , <J?b and V. Again, we refer reader to AppendixEIfor our notation. We treat their 
masses as free parameters and investigate the possibility for consistent scenario with the exact 
one-loop unification. Since all other fields in the Higgs sector, i.e., ^t, S§ and <£ c , simply worsen 
unification we simply assume they live at or above the grandunifying scale. 

In order to present consistent analysis we now discuss the constraints coming from proton decay 
on Bij coefficients. These enter via Eq. (|5b[) and assumption that My = Mqut- As we show these 
constraints are rather weak if the gauge d = 6 contributions are dominant as is usually assumed 
in non-supersymmetric GUTs Q]. For example, if we use the latest bounds on nucleon decay 
lifetimes we obtain, in the context of an SU (5 non-supersymmetric GUT, in the case of maximal 
(minimal) suppression in the Yukawa sector |l9| My > 2.5 x 10 13 GeV (My > 1.5 x 10 15 GeV). (The 
minimal suppression case corresponds to the GG scenario with Yjj = Y^ and Yd = Yg, where Yu, 
Yd and Ye are the Yukawa matrices of charged fermions. Non-renormalizable contributions violate 
both of those relations. The same is also true for the running in the Yukawa sector from the GUT 
scale where those relations hold to the scale relevant for the Yukawa couplings entering nucleon 
decay. On the other hand, maximal suppression corresponds to a case with particular relation 
between unitary matrices responsible for bi-unitary transformations in the Yukawa sector 19] that 
define physical basis for quarks and leptons.) In both cases we use cxq UT = 35 and the best limit 
on partial lifetime which is established for p — > 7r°e + decay channel (r > 5.0 x 10 33 years). This 
gives conservative bound for the suppressed case since it is always possible to rotate away proton 
decay contributions for individual channels [lfl| . 

The uncertainty in extracting the limits on My from experimental data is easy to understand. 
Namely, even though the nucleon lifetime is proportional to My, which would make extraction 
rather accurate and precise, the lifetime is also proportional to the fourth power of a term which 
is basically a sum of entries of unitary matrices which are a priori unknown unless the Yukawa 
sector of the GUT theory is specified and which, in magnitude (see Eqs. (|B2[H. can basically vary 
from Vub to 1 18, 3]. (For full discussion see also Appendix^!) If we now adopt the Mqut = My 
assumption and use Eq. (|5b|) the above limits translate into B12 < 7.0 (B12 < 6.1) for the suppressed 
(unsuppressed) case. So, all SU (5) GUTs with B12 > 7.0 are excluded by the usual gauge d = 6 
contributions to proton decay. The theories with 6.1 < By± < 7.0 require "special" structure of 
the Yukawa sector; the closer the B12 to the upper limit is the more "special" structure is needed. 
Finally, any SU(5) GUT with B12 < 6.1 has not yet been probed by proton decay experiments. 
(Again, this is all based on the one-loop analysis. Any more accurate and precise statement must 
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be based on the two-loop treatment with a proper inclusion of the threshold corrections.) 

In order to avoid problems with proton decay without requiring too much conspiracy in the 
Yukawa sector we pursue the solutions where the superheavy gauge bosons are as heavy as possible. 
So, how heavy can they be given the particle content of the SU(5) scenario with an extra 15 of 
Higgs? In order to answer that we first naively set masses of S3, <E» a and to Mz- This in turn 
yields the lowest possible value of B12 to be 6.4 (6.33) for n = 1 (n = 2) which translates via 
Eq. (HQ) into M GUT = 3.2 x 10 14 GeV (M GUT = 4.4 x 10 14 GeV). (We include n = 2 case in our 
considerations since it might be relevant in addressing Baryon asymmetry of the Universe.) In 
this naive estimate My is either equal to or slightly above Mqut- From the previous discussion 
on the My limits we see that there is a need for small suppression in Yukawa sector in order to 
satisfy experimental limits on proton decay. This suppression, as we show later, amounts to a 1/5 
of the unsuppressed case. When compared to the available suppression (~ 1/V u b) it comes out to 
be around 2%, which can easily be accomplished. 

Note, however, that in order to have exact unification crucial thing is to satisfy Eq. (|5a[). Thus, 
it is better to ask for which mass spectrum that satisfies Eq. (|5a|) we obtain the highest possible 
value for My or equivalently the smallest possible value for Byi- As it turns out the answer to 
this question is unique within our scenario. To show that we first assume that the relevant degrees 
of freedom that improve the running, i.e., S3, <3? a and <$>b, contribute in pairs, e.g., a degenerate 
pair (S3,<I> a ) is light and is at Mqut, and treat only the n = 2 case. With those constraints 
we generate three possible combinations which yield results summarized in Table HT1 (We address 
both the n = 1 and n = 2 case at the two- loop level in Appendix O) 

TABLE II: AB23 and ABi 2 corrections due to the degenerate pairs of fields and the associated scales for 
the n = 2 case. 

(S 3 ,$ a ) (S 3 , $6) (*a,$b) 

AB 23 fr + f |r + | |r + | 



A^i 2 -&r-& -Afr-A 

M GUT 5 x 10 13 GeV 9 x 10 13 GeV 



M r 1 TeV 200 GeV 



• The (S3, $ a ) case with n = 2 exactly mimics the n = 8 case in terms of quantum numbers. 
(Recall that it takes at least eight light Higgs doublets on top of the Higssless SM content to 
unify the couplings. Associated corrections to the Higssless SM coefficients are AB23 = | r +| 
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and AB12 = —j^f — j*;, where r, as denned in Eq. ©, is very close to one and we take 
two of the doublets to be at Mz-) The unification scale is rather low and very close to the 
experimentally set limit for maximally suppressed case. Lightness of <3? a goes against the idea 
behind the type II see-saw if one assumes that the parameter C3 (see Eqs. (jA5|) and (|AlljO 
is at the GUT scale, but at this point the scenario is not ruled out experimentally. 

• The ($ a , case has a slightly higher unification scale than the (X^,^) case. This time 
both <I> a and have mass in phenomenologically interesting region. Lightness of $ a again 
requires large suppression in the Yukawa sector for neutrinos to generate correct mass scale 
via type II see-saw. However, such a suppression would be beneficial in suppressing novel 
contributions to proton decay due to the mixing between and ^t- 

• The (E3, <!>;,) case is the most promising. Even though it fails to unify at the one-loop 
level its correction to B12 is the largest of all three cases. As such, it represents the best 
possible candidate to maximize My. Moreover, the $ a contribution to the running to produce 
unification for light S3 and &b is small which implies that its mass could be in the range 
that is optimal for the type II see-saw for the most "natural" value of C3 coefficient. Again, 
the (£3,<&&) case not only maximizes My but also places M$ a at the right scale to explain 
neutrino masses. 

The three special cases discussed above all demonstrate that large My scale prefers §b light 
regardless of the relevant scale of other particles since it is <3?6 coefficients that decrease B12 the 
most. These conclusions persist in more detailed one- and two-loop studies. 

Is there a way to tell between the three limiting cases we just discussed? The (£ 3 ,<I> a ) case 
can be tested and excluded by slight improvement in the nucleoli lifetime data; other low energy 
signatures depend on how light <I> a is. One could also test and distinguish between the (23,$^) 
and ($ a , cases since both favor light &b leptoquarks that can be detected by LHC. If and when 
these are detected the two cases could be distinguished by the scalar leptoquark contributions 
to the rare processes. In the ($ a ,<^b) case the suppression in the neutrino Yukawa sector would 
selectively erase some of these contribution while in the (£3, case all these contributions would 
be sizable. Expected improvements in the table-top experiments would then be sufficient to tell 
the two. 

Since we have Mqut and masses of V, S3, $ a and <3?b as free parameters and only two 
equations — Eqs. (|5a[) and l)5bjl — we present four special cases based on certain simplifying as- 
sumptions in Figs, n ID and El and discuss each case in turn. All examples we present generate 
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consistent unification in agreement with low energy data. (Note that the change in the parameters 
also affects the value of ot^fj T . We do not present that change explicitly, which, for the range of 
values we use, vary from 36 to 40. In our plots we also allow My to be at most factor of three 
or four lower than the GUT scale. Once we switch to two-loop analysis with threshold corrections 
accounted for we appropriately set My = Mqut-) 



23456789 




23456789 
M GUT (10 14 GeV) 



FIG. 1: Plot of lines of constant value of My (solid lines) and Af$ Q (dashed lines) in the Mgut~{Ms 3 = M$ b ) 
plane. The area to the left of a steep solid line denotes the region where My does not contribute to the 
running, i.e., My > Mgut- To generate the plot we consider exact one-loop unification and use central 
values for the gauge couplings as given in the text. This is a scenario with one light Higgs doublet (n = 1). 

In Fig. ^ we present the n = 1 case when the pair (£3,$;,) is taken to be degenerate with the 
mass close to electroweak scale (~ O(10 2 ) GeV). The parameter C3 has to be between 10 3 and 10 7 
GeV to explain the neutrino mass through the type II see-saw if the Yukawa coupling for neutrinos 
is of order one. On the other hand, the gauge boson mass varies only slightly for a given range of 
M$ b and Ms 3 around 2.5 x 10 14 GeV. Clearly, unification itself allows M$ 6 and M^ 3 to be much 
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heavier then 1 TeV on account of decrease of M$ but in that case My would be getting lighter. 
This, on the other hand would require additional conspiracy in the Yukawa sector in order to 
sufficiently suppress proton decay to avoid the experimental limit. 

The two light Higgs doublet case is presented in Fig. [5J This case is well motivated on the 
Baryogenesis grounds. Namely, the interaction of the 15 of Higgs explicitly breaks B — L symmetry 
(see Appendix^). This opens a door for possible explanation of the Baryon asymmetry in the 
Universe within our framework. However, since the successful generation of Baryon asymmetry 
requires at least two Higgses in the fundamental representation we study a consistent unification 
picture for that case. (See references for the Baryogenesis mechanisms in the context of SU(5) 
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22, 231). The n = 2 case has 



model with two Higgses in the fundamental representation 
higher scale of superheavy gauge bosons compared to the n = 1 case (this does not hold at the 
two-loop level though) and the mass of the field $ a is in the region relevant for the type II see-saw. 
Hence, if the mass of scalar leptoquarks is in phenomenologically interesting region (~ O(10 2 ) GeV) 
we can explain neutrino masses naturally. Again, as in the n = 1 case, the gauge boson mass varies 
very slightly, this time around 3 x 10 14 GeV. Given the last two examples we can again conclude 
that the exact unification in this minimal realistic scenario points towards light scalar leptoquarks. 

In order to understand better these results we show two more examples in Fig. |3J This time we 
set S3 = <5 a for simplicity and present both n = 1 and n = 2 cases. This scenario is disfavored 
by the fact that My tends to be "small" but cannot be excluded on experimental grounds. It is 
evident from the plot that My does not depend on a number of light Higgs doublets and Mqut- 
The reason for that is very simple and is valid only at the one- loop level. Namely, the ratio B23/B12 
is the same for typ coefficients (= —5/2) as for the sum of corresponding S3 and <£ a coefficients 
provided these are degenerate. Thus, any change in the number of light doublets in Eq. ()5a|) is 
simply compensated by the change in degenerate mass of S3 and $ a fields for a fixed value of <&b 
mass. This trend can be clearly seen in Fig. El The mass of $ a is generally rather low to generate 
neutrino mass of correct magnitude unless Yukawa couplings of neutrinos are extremely small. 

In all our examples, S3 is allowed to be much lighter than Ss. But, this seems in conflict with 
the treedevel analysis of the S potential that is invariant under S —* — S transformation which 
yields a well known relation M^ 3 = 4Ms g |l7] • This apparent mismatch has a simple remedy. 

In order to generate sufficiently large corrections to charged fermion masses via higher- 
dimensional operators in the Yukawa sector we need terms linear in S/Mp;. If that is the case it 
is no longer possible to require that the Lagrangian is invariant under transformation S — > — S. 
It is then necessary to include a cubic term into the S potential besides the usual quadratic and 
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FIG. 2: Plot of lines of constant value of My (solid lines) and Af$ Q (dashed lines) in the Mqut-{Ms 3 = M$ b ) 
plane. The area to the left of a steep solid line denotes the region where My does not contribute to the 
running, i.e., My > Mqut- This is the n = 2 scenario. 

quartic ones. But, the potential with the cubic term (TrS 3 ) violates the validity of Ms 3 = 4Ms 8 



relation 
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261 ] and allows a possibility where S3 is light while Ss is superheavy. We analyze 
this situation in Appendix^] in some detail. Note that we do not require nor insist on the lightness 
of S3 though. In Appendixinjwe present the two-loop analysis of the scenario where <3? a is relatively 
light (~ 10 7 GeV) and S3 is at the GUT scale. Our intention is solely to demonstrate that there 
are more possibilities available unless additional assumptions, such as S — > — S transformation, are 
imposed on the SU(5) theory. Note, however, that the maximization of My always requires <]?{, to 
be very light (~ 10 2 GeV). In Fig. 0] we show an example where it is possible to achieve unification 
at the two-loop level (See details in AppendixEJ) for n = 1, M$ b = 250 GeV, M$ a = 1.54 TeV and 
the field S3 is at the GUT scale. 

What about the possible mass spectrum of <£ a , and $ c ? The relevant potential is in Ap- 
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FIG. 3: Plot of lines of constant value of My (solid lines) and M$ = Ms 3 in GeV units (dashed lines) in 
the McuT-(M$ b ) plane. The scenario with one (two) Higgs doublet(s) corresponds to long (short) dashed 
lines. 

pendix^J Clearly, there are more parameters than mass eigenvalues. The tree- level analysis revels 
that it is possible to obtain any possible arrangement including, for example, << & a < <£ c . 
This sort of split is quite similar to the split behind the well-known Doublet- Triplet problem. 

Our framework yields rather low mass for vector leptoquarks that varies within very narrow 
range around 3 x 10 14 GeV for a most plausible scenarios. This makes the framework testable 
through nucleon decay measurements. (More precisely, large portion of the parameter space of 
the setup has already been excluded by existing measurements on nucleon lifetime.) It is easy 
to understand this generic and robust prediction. We have seen that Mqut can be at most 
3.2 x 10 14 GeV (4.4 x 10 14 GeV) if we exclude the V contribution from the running for the n = 1 
(n = 2) case. If we now start lowering My below Mqut we lower B\2 as well. This, on the 
other hand, starts to increase Mqut but still keeps My at almost the same value. Basically, the 
"decoupling" of My and Mqut takes place, where the mass of X and Y gauge bosons remains 
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FIG. 4: Unification of the gauge couplings at the two-loop level for central values of low-energy ob- 
scrvablcs 15]. The SM case with n — 1 is presented by dashed lines. Solid lines correspond to the n = 1 
scenario with <&& and <I> Q below the GUT scale. Vertical lines mark the relevant scales: Mz, M$ b = 250 GeV, 
AU a = 1.54 TeV and M GUT = 0.96 x 10 14 GeV. 

in vicinity of the value before the "decoupling" while Mqut rapidly approaches the Planck scale. 
Another way to say this is that the Bij coefficients of superheavy gauge fields are very large 
compared to all other relevant coefficients (see TableCJ). Thus, any small change in My corresponds 
to a large change in other running parameters. 

Let us finally investigate possible experimental signatures coming from our consistent minimal 
realistic SU (5) model. 

IV. EXPERIMENTAL SIGNATURES 

Our framework has potential to be tested through the detection — direct or indirect — of light 
leptoquarks and/or observation of proton decay. Let us investigate each of these tests in turn. 

• Light leptoquarks 

To get consistent unification in agreement with low energy data, neutrino mass and proton 
decay in our minimal framework we generate very light leptoquarks <£>{,. The lighter the 
is the heavier the V becomes. Thus, in the most optimistic scenario M$ b is close to the 
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present experimental limit ~ O(10 2 ) GeV. In what follows we specify all relevant properties 
of <&b and existing constraints on its couplings and mass. 

The 5^C/i a fe$5;, coupling yields the following interactions: 

<f T Ch§ b l = d cT Ch{4>\e- $v), (6) 

where the leptoquarks (f>\ and have electric charges 2/3 and —1/3, respectively, and 
symmetric matrix h coincides with the Yukawa coupling matrix of Majorana neutrinos (= Y v ) 
if we neglect the Planck suppressed operators. (See the last line in Eq. (|A4[) .) The above 
leptoquark interactions in the physical basis read as: 

<f T C (D^E*)K*(Vf MNS rY^(V^ MNS YK* 4 e, (7) 



d cT C {D T C E*)K*(V^ MNS YY^ <pl „, (8) 

where Dq and E are the matrices which act on d quarks and e, respectively, to bring them 
into physical basis. (See Appendix El for exact convention.) K is a matrix containing three 
CP violating phases, and Vp { MNS is the leptonic mixing in the Majorana case. (In the GG 
SU (5) where Ye = Y^ one has Dc = E. However, that is not the case in a realistic model 
for fermion masses.) 

There are many studies about the contributions of scalar and vector leptoquarks in different 
For a model independent constraints on leptoquarks from rare processes 



processes |27j, |2j 
see for example 



28, 



29l | . The most stringent bound on the scalar leptoquark coupling to 
matter comes from the limits on fi-e conversion on nuclei j^] . The bound we present should 
be multiplied by (M$ b /100 GeV) 2 . In our case it reads (D^hE)n(E^ D^) 2 i < 1(T 6 . The 
bounds for all other elements of (DghE)ij(E* h/ Dg)ki and (D^hN)ij(N*h^ Dq)^ are weaker. 

The currents bounds on leptoquarks production are set by Tevatron, LEP and HERA 
Tevatron experiments have set limits on scalars leptoquarks with couplings to eq of Mlq > 
242 GeV. The LEP and HERA experiments have set limits which are model dependent. 
The search for these novel particles will be continued soon at the CERN LHC. Preliminary 
studies by the LHC experiments ATLAS 32] and CMS [3^] indicate that clear signals can 
be established for masses up to about Mlq ~ 1.3 TeV. For several studies about production 
of scalar leptoquarks at the LHC, see Ref. [3^ . Thus, it could be possible to test our 
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scenario at the next generation of colliders, particularly in the Large Hadron Collider (LHC) 
at CERN, through the production of light leptoquarks. Therefore even without the proton 
decay experiments we could have tests of this non-supersymmetric GUT scenario. 

Proton decay 

Proton decay is the most generic prediction coming from matter unification; therefore, it 
is the most promising test for any grand unified theory. (For new experimental bounds 



sec 



35j.) In our minimal and consistent scenario the relevant scale for gauge bosons is 
around 3 x 10 14 GeV regardless of how high the GUT scale goes in order to get consistent 
unification in agreement with low energy data. Careful study within the two-loop context 
with the inclusion of threshold effects revels that the highest possible value of My in the 
n = 2 (n = 1) case is 3.19 x 10 14 GeV (3.28 x 10 14 GeV) for central values of coupling 
constants 15] while the 1 a departure allows for the maximum value of 3.35 x 10 14 GeV in 
the n = 2 case, for example. 

There are several contributions to the decay of the proton in our minimal scenario. We have 
the usual Higgs and gauge d = 6 operators but there are also new contributions due to the 
mixing between ^>t an d with <& b being extremely light in our case. These contributions 
are very important. Using the relevant triplet interactions: 

q a C A ab q b m T + qa C C ab l b ^ + u c a C D ab <f ^ + u c a C B ab e% tf T) (9) 

(for the expressions of A, B_, C_, and D. matrices see for example ,11]) and the interaction 
term 5^ Ch ab <&5 b , it is easy to write down the contributions for the B — L non-conserving 
decays p — > {K + , 7r + , p + )vi, and n — > (jr°, p°, rf , u>°, K )^. We present the relevant diagram 
in Fig. El Q- 
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FIG. 5: Contributions to the decay of the proton induced by the 15 of Higgs. 

Notice that in this scenario we have the usual B — L conserving decays, i.e., the decays into 
a meson and antileptons, and the B — L non-conserving decays mentioned above. Since 
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has to be light, the B — L violating decays are very important. The rate for the decays into 
neutrinos is given by: 

r(p -► K + v) = Vrfp^ K + Vi) = A^A\ 3 . . x 

4m 2 D 2 

\pC(vi,s,d c ) + dCKs c /)| 2 -V 

+ \pC^,d,s G ) + aCkrf c lS c )| 2 [l + mp( f m + 3F) ] 2 (10) 

where C(v i: d a ,d c ) = (U T {A + A T )D) la and C(ui,d G ,d c ) = (D^D^U^ai- (See Appendix 
B for notation.) 

As you can appreciate from the above expressions, the predictions coming from these con- 
tributions are quite model dependent. Using m p = 938.3 MeV, D = 0.81, F = 0.44, 
m B = 1150 MeV, f n = 139 MeV, A L = 1.43, and a = (3 = a = 0.003 GeV 3 we get: 

M 4 M 4 

1.95 x 10" 64 GeV" 6 % * b > [0.19 \C(^,s,d c ) + C(u h s c , d c )\ 2 + 

4 

+ 2A9\C(^,d,s c ) + C(vi,d c ,s c )\ 2 ] (11) 

Let us see an example, using the values My T = c 3 = 10 14 GeV and M$ 6 = 10 3 GeV, the 
left-hand side of the above equation is equal to 1.95 x 10~ 24 ; therefore, the sum of the C 
coefficients has to be basically 10 -12 . In the case that coefficient C3 is smaller (~ 10 6 GeV), 
a possibility that is not excluded, the sum of C coefficients would be around 10 -4 , which 
is their "natural" value. Moreover, the scenario would then prefer $ a at the same scale 
(~ 10 6 GeV) if Y v is taken to be proportional to Y e . Also we can suppress the relevant 
contributions in different ways. For example, we could choose = —A^ and = 
except for i = j = 3, or set to zero these coefficients in specific models for fermion masses. 

In any case, if the gauge d = 6 contributions are the dominant ones for proton decay, we can 
get the following bounds for the proton lifetime (see Appendix [BJ) allowing the full freedom 
in the Yukawa sector: 

1 x 10 31 years < r p < 2 x 10 38 years. (12) 

Here we use ot^j T = 39 and M GUT = 3.2 x 10 14 GeV. (See Fig. © for example.) Having 
in mind the experimental limit of 5.0 x 10 33 years |l5l | we see that a significant portion of 
available parameter space has already been excluded. We hope that in the next generation 
of proton decay experiments this scenario would be constrained even further. 
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How does this scenario compare to other possible extensions of the GG model? We mention 
only few listing them by increasing order in particle number. 

• The most obvious extension is to add one more fundamental representation in the Higgs 
sector to the n = 2 scenario we analyze the most at the one-loop level. This addition would 
not raise but actually lower My since the scalar leptoquarks which influence Byi the most 
would get slightly heavier than in the n = 2 case. In certain way, this actually makes the 
n = 1 scenario with the 15 of Higgs very unique. It is the minimal extension of the GG 
model with the highest available scale for My. We focus our attention on the n = 2 case 
on the grounds of Baryogenesis. In the same manner, the n = 3 scenario would be well 
motivated by the possibility of addressing the issue of the SM model CP violation |3J3] (see 
also |3] and references therein). 

• Very interesting possibility would be the n = 1 case with two 10s of Higgs. Such a scenario 
could have a very high GUT scale and still very promising phenomenological consequences 
due to light leptoquarks. For example, successful two-loop unification with light ^s (there 
are two now) at 250 GeV requires M S3 = 2.1 x 10 11 GeV and the GUT scale at 1.0 x 10 15 GeV 
for central values of ci^s at Mz- This model would also require at least right-handed neutrinos 
and non-renormalizable operators to be completely realistic. 

• Another possibility would be the n = 1 case with one 10 and one 15 of Higgs. Such a 
scenario would have a same maximal value for the GUT scale as the two 10 case above. To 
be completely realistic it would require non-renormalizable operators. 

• The next scenario is the one proposed by Murayama and Yanagida [la] (MY). They shown 
that addition of two 10s of Higgs to the GG model and with n = 2 it is possible to achieve 
unification for extremely light scalar leptoquarks in the 10s. This allows MY to forward 
a "desert" hypothesis, within which the particles are either light (~ O(10 2 ) GeV) or heavy 
(~ O(10 14 ) GeV). (Note that the scalar leptoquarks with exactly the same quantum numbers 
as $b (c 15) reside in the 10 as well. To see that one can use 5® 5 = 10© 15. Note that the 
couplings of $b in the 10 (15) to the 5s are antisymmetric (symmetric).) Their model is ruled 
out by direct searches for scalar leptoquarks due to the extreme lightness of <3?ftS. However, 
if one allows for splitting between £3 and Ss one can raise unification scale sufficiently to 
avoid proton decay bounds and resurrect their model although one would have to abandon 
the "desert" hypothesis MY forwarded. The additional contribution to the Higgsless SM 
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coefficient in this case is basically AB12 = —jgr$ b — y^E 3 — j§ which is much more than 
any of the cases considered in Table [0] For example, for M§ b = 300 GeV we find, at the 



one-loop level, M Sa = 10 10 GeV and M GUT = 9 x 10 14 GeV. Note that the presence of the 
10 of Higgs yields the same type of coupling as we specify in Eq. (J5J except that h, in this 
case, would be antisymmetric. This, however, would not a priori prevent proton decay. If 
there are two or more Higgses in the fundamental representation in the model there exist 
the proton decay process schematically represented in Fig. |1] unless additional symmetry 
is introduced. Finally, this model would also require at least right-handed neutrinos and 
non-renormalizable operators to be completely realistic. 

The Georgi-Jarlskog model [lfl| with an extra 45 of Higgs to fix fermion masses represents 
natural extension of the GG model. However, even though unification takes place ^ ll^t 
and the extension has good motivation the predictivity of the model is lost unless additional 
assumptions are introduced. More minimal extension than that would be, for example, an 
extra 24 of Higgs to the scenario we consider. 



V. SUMMARY 



We have investigated the possibility to get a consistent unification picture in agreement with 
low energy data, neutrino mass and proton decay in the context of the minimal realistic non- 
supersymmetric SU(5) scenario. This scenario is the Georgi-Glashow model extended by an extra 
15 of Higgs. As generic predictions from the running of the gauge couplings we have that a 
set of scalar leptoquarks is light, with their mass, in the most optimistic case, being around 
O(10 2 -10 3 ) GeV. This makes possible the tests of this scenario at the next generation of collider 
experiments, particularly in the Large Hadron Collider (LHC) at CERN. In the "least" optimistic 
scenario the mass of the scalar leptoquarks would be around 10 5 GeV. The proton decay issue 
has been studied in detail, showing that it is possible to satisfy all experimental bounds with 
very small, at 2% level, suppression in Yukawa sector. Rather low scale of vector leptoquarks 
(~ 3x 10 14 GeV) already allows for significant exclusion of available parameter space of our scenario. 
Further reduction is expected in near future with new limits on proton decay lifetime. We have 
particularly studied the case with two Higgses in the fundamental representation since in this case 
it could be possible to explain the Baryon asymmetry of the Universe. We have also compared this 
scenario with other, well motivated, extensions of the GG model. There are uncertainties related 
to predictions of the proposed scenario but, in view of the fact that it truly represents the minimal 
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realistic extension of the GG model, the same is even more true of all other extensions. 
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APPENDIX A: PARTICLE CONTENT AND RELEVANT INTERACTIONS 

In this appendix we define a minimal realistic non-supersymmetric SU(5) model. By min- 
imal we refer to the minimal number of physical fields such a model requires. In the GG 
model the matter is unified in two representations: 5 a = (d c ',l) a , and 10 a = (u c , q, e c ) a , 
where a = 1,2,3 is a generation index. The Higgs sector comprises 5h = = (^d^t) an d 
24h = S = (Sg, S3, S( 3) 2), S(3 j2 ), S24). However, this model does not achieve unification and fails 
to correctly accommodate fermion masses; therefore, it is ruled out. In the introduction we discuss 
how it is possible to solve all phenomenological problems of the GG model introducing a minimal 
set of fields. Namely, it is sufficient to introduce an extra representation: 15h = $ = (<I> a , $5, <I> C ). 
The SM decomposition of the Higgs sector is given by: 

24 H = £ = (8, 1, 0) + (1, 3, 0) + (3, 2, -5/6) + (3, 2, 5/6) + (1, 1, 0), (Al) 
15 H = $ = (1, 3, 1) + (3, 2, 1/6) + (6, 1, -2/3), (A2) 
5 H = * = (1, 2, 1/2) + (3, 1, -1/3). (A3) 

The relevant Yukawa potential, up to order 1/Mpi, is 

io^io^" 1 + loiViabiof + 10^f 2ab 10 k b n ¥j^-J 

+ v>W h + i^Mggi , n (A4) 

Mpi 

where i,j,k,l,m represent SU(p) indices. Impact of the non-renormalizable operators on the 
fermion masses is discussed in 9]. Notice that we could replace Mpi by a scale A, where Mqut < 
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A < Mpi, if we do not assume a desert between the GUT scale and the Planck scale. (See for 
example reference 3911 .1 

nn 

The Higgs scalar potential, manifestly invariant under SU(5), is (see for example |4fll I41l|): 
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(A5) 



We assume no additional global symmetries. It is easy to generalize this potential to describe 
the case of two or more Higgses in the fundamental representation. (We do not include non- 
renormalizable terms in the Higgs potential since the split between E3 and Eg masses that we 
frequently use in running can already be achieved at the renormalizable level.) 

The condition that the symmetr y b reaking to the SM is a local minimum of the Higgs potential 



for S (the first line in Eq. ()A5 
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where dimensionless variables are defined as 
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The vacuum expectation value of E is (E) = A/V30 diag(2, 2, 2, —3, —3), where 
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Finally, the mass of X and Y gauge bosons is given by 



My = \j — gcurX 



(A7) 
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Here we note the following. In the limit that cs -> we obtain the well know results: M2 



4M| g Q. However, in the limit where A — > f 2 - 1 ^ we obtain 



£3 



V vr o-gut J 



Ml 3 -> 0, 



7 "5 
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Clearly, it is technically possible to achieve a large split between Ms 3 and Ms 24 although this is 
highly unnatural. 

The relevant interactions for the see-saw mechanism are the following: 

^see-saw = -AfJ, Tr <5>l$ a -Y v l T C <S> a I + c 3 + h.c. (A9) 

where: 

*a=^_^^j, i=r|, *S = (fi°,-a r+ )- (Aio) 

In this case the neutrino mass is given by: 



M v « n c 3 -g (All) 



<J/°> 2 

Ml 

<Pa 

This is | he so-called type II sec-saw I. 5|. Notice that from this equation in order to satisfy the 
neutrino mass experimental constraint M^Jo$ has to be around 10 13-14 GeV. 

APPENDIX B: PROTON LIFETIME BOUNDS 

The dominant contribution towards nucleon decay in non-supersymmetric GUT usually comes 
from the gauge d = 6 proton decay operators. Even though these operators carry certain model 
dependence we have recently shown that unlike in the case of d = 5 operators we can still establish 
very firm absolute bounds on their strength that are equally valid for all unifying gauge groups . 
We investigate the impact these bounds have on the non-supersymmetric grand unified model 
building in what follows. 

The upper bound for the total nucleon lifetime ^| in grandunifying theories, in the Majorana 
neutrino case, reads 

M 4 

t p < 1.1 x 10 41 years^-^, (Bl) 
a GUT 
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where My — the mass of the superheavy gauge bosons — is given in units of 10 16 GeV. We stress 
that there exist no upper bound for partial lifetimes since we can always set to zero the decay rate 
for a given channel. 

In order to fully understand the implications of the experimental results we now specify the 
theoretical lower bounds on the nucleon decay in GUTs. These are applicable in the case of non- 
supersymmetric models if the gauge d = 6 contributions are the dominant as well as in the case 
of supersymmetric models where the d = 4 and d = 5 operators are either forbidden or highly 
suppressed. 

The relevant coefficients for the proton decay amplitudes in the physical basis of matter fields 
are Q: 

c(e°,dp) = k\ [ V?Vf + {ViV UD )^(V 2 V^ D r\ (B2a) 
c(e a ,d%) = fcf V?Vf+ k 2 2 (ViV^fHViVuDVlVz) 1 ", (B2b) 
c(v h d a , d c p ) = k\ {V x V UD ) la {V z V EN f l + k\ Vf a {V 1 V UD vlv s V EN ) n , a = 1 or = ^B2c) 
c(vf,d a ,dC) = k\ \(y A vb D f\ul N V 2 ) la + vf a (U EN V 2 v£ D ) 11 }, a = 1 or (5 = 1, (B2d) 

where a, (3 = 1,2. The physical origin of the relevant terms is as follows. The terms proportional 
to ki(= gcuTMy 1 ) are associated with the mediation of the superheavy gauge fields V = (X, Y) = 
(3, 2, 5/6), where the X and Y fields have electric charges 4/3 and 1/3, respectively. This is the case 
of theories based on the SU(5) gauge group. On the other hand, an exchange of V' = (X',Y') = 
(3, 2, —1/6) bosons yields the terms proportional to k 2 (= gcuTMy}). In SO(W) theories all these 
superheavy fields are present. 

The relevant mixing matrices are Vi = U j c U, V 2 = E j c D, V 3 = D j c E, V 4 = D j c D, V UD = WD, 
V E n = E^N, U E n = E ci N c , and V UD = U^D = KiV C kmK 2 , where K x and K 2 are diagonal 
matrices containing three and two phases, respectively. The leptonic mixing Ven = ^^pmnS^ 
in case of Dirac neutrino, or Ven = K-zVpMNS in * ne Majorana case. Vp MNS and Vp M MNS are 
the leptonic mixing matrices at low energy in the Dirac and Majorana case, respectively. Our 
convention for the diagonalization of the up, down and charged lepton Yukawa matrices is specified 
by U^YuU = Yff g , D T C Y D D = Y* ag , and E%Y E E = Y E ^. 

To establish the lower bound on the nucleon lifetime we first specify the maximum value for 
all the coefficients listed above for 50(10) theory only. The SU(5) case is well known and can be 
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reproduced by setting k,2 = in the expressions below. The upper bounds are 

<<%,dp)so{iO) < 2 hf, (B3) 
50(10) < k\ + kl (B4) 

3 

Y J c^i,d a ,d^y som c(iy l ,d 1 ,df) soiW ) < kf5 ps + kl + 2klkl (B5) 
i=i 

3 

^c(vf,d a ,d%y som c(vf,d 7 ,d$) som < kl[3 + 5 a % (B6) 

i=l 

which translates into the following bounds on the amplitudes in the case the neutrinos are Majorana 
particles: 

r(P - - + -) < Sl^il«| 2 (l + O + Ffikf + A|] 2 , (B7) 



Svr/ 2 ' 

(m 2 p -m l K f „ 2 rr 2m p D 1 



r ^* +p > £ ^7>|-^ a i 3x f[iS : i 2+ [ 1 + S (D+3F)l2lx (B8) 

(*? + fcl) 2 + + + aF)][2*?*§ + fcf}, (B9) 

r(p^7r°e+) < -^Al\a\ 2 (l + D + F) 2 [hk A 1 + 2klkl + kil (BIO) 



r(p - K%) < " p 'f A\\a\ 2 [l + ^(D - F)] 2 [5fc 2 + 2k\k\ + fc 2 4 ]. (Bll) 



{ m 2 p -m 2 K ) 2 2 2 m f 
87rf4m^ m B 

Using these expressions it is easy to extract lower bounds on lifetimes. In Table ITTTl we list all lower 
bounds for the proton lifetime in SU(5) and SO(10) models. Again, we use m p = 938.3 MeV, 
D = 0.81, F = 0.44, m B = 1150 MeV, f n = 139 MeV, A L = 1.43, and the most conservative value 
a = 0.003 GeV 3 . In the case of 5O(10) models we set k% = &2 for simplicity. Note that lower 

TABLE III: Lower bounds for partial proton lifetime in years for the Majorana neutrino case in units of 
My/aQ UT , where the mass of gauge bosons is taken to be 10 16 GeV. 

Channel T J f C/( ' 5 ^ Tp 01 " 10 ^ 

P^tt+P 7.3 x 10 33 1.8 xlO 33 

p -> K+v 17 .4 x 10 33 4.8 x 10 33 

p^7r°e+ 3.0 xlO 33 1.8 xlO 33 

p -> K°e+ 8.5 x 10 33 5.3 x 10 33 



bounds are well defined for the partial lifetimes while upper bound is meaningful for the total 
lifetime only. 

Finally, we can establish the theoretical bounds for the lifetime of the proton in any given GUT. 
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In what follows we use 



10 34 years < t p < 10 41 years 



a GUT a GUT 



(B12) 



These bounds are useful since we can say something more specific about the allowed values of My 
or olqut or both. For example, if we take acuT we can put lower limit on the value of My using 
experimental data on nucleon lifetime. Also, given the value of My and olqut we can set the limits 
on the proton lifetime range within the given scenario. 

APPENDIX C: THE TWO-LOOP RUNNING 



We present the details on the two-loop running. In order to maximize My one needs extremely 
light <]?{,. In any given mass splitting scheme we thus set the mass of at 250 GeV which is 
just above the present experimental limit of 242 GeV. Then, S3 or $ a or both are allowed to vary 
between M$ 6 and Mqut in order to yield unification and all other fields except for the SM ones are 
taken to be at the GUT scale. The identification My = Mqut is justified through the inclusion of 



boundary conditions at the GUT scale 



43] 



«i \GUT = a GUT 



12vr : 



{Ai,A 2 ,A 3 } = {5,3,2}, 



(CI) 



and the relevant two-loop equations for the running of the gauge couplings take the well-known 
form 

Q 

dai(/j,) h 



d /x 



2vr 



(C2) 



bi and bij coefficients for the SM case for arbitrary n are well-documented (see |44j for general 
formula). In addition to those we have: 
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which should be added to the SM ones at the appropriate particle mass scale. The two-loop 
Yukawa coupling contribution to the running of the gauge couplings (with the corresponding one- 
loop running of Yukawas) is not included in order to make meaningful comparison between n = 1 
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and n = 2 cases. (For the n = 2 case one vacuum expectation value of the light Higgs doublets is 
arbitrary and needs to be specified in order to extract fermion Yukawas for the running, i.e. tan (5 
ambiguity. There is no such ambiguity present in the n = 1 case since the only low-energy vacuum 
expectation value is accurately determined by electroweak precision measurements.) 
The outcome of the exact numerical unification is presented in Figs. El [Zl and [HI • 
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FIG. 6: Unification of the gauge couplings at the two-loop level for central values of low-energy ob- 
servables The SM case with n = 2 is presented by dashed lines. Solid lines correspond to the n = 2 
scenario with <!>{, and S3 below the GUT scale. Vertical lines mark the relevant scales: Mz, M$ b = 250 GeV, 
M Ss = 4.95 x 10 4 GeV and M GUT = 3.19 x 10 14 GeV. 

The GUT scale as well as appropriate intermediate scales are indicated on the plots. For 
example, in the n = 2 scenario with light S3 in Fig. the GUT scale is close to the one-loop 
results (see Fig. particular) and comes out to be 3.19 x 10 14 GeV for central values of coupling 
constants jis| . The la departure allows for the maximum value of 3.35 x 10 14 GeV in that case. 
Note that the n = 1 case with light S3 presented in Fig. [7| yields somewhat higher GUT scale. 
The reason behind this trend is simple: the §b contribution to Byi is seven times that of the Higgs 
doublet but, at the same time, they contribute the same to I?23- Thus, is more efficient in 
simultaneously improving the running and raising the GUT scale than any extra Higgs doublets. 

One might object that lightness of S3, which requires miraculous fine-tuning, makes the scenario 
with an extra 15 rather unattractive. Again, we do not insist on S3 being at intermediate scale. 
Note that the unification with almost the same GUT scale as in the intermediate S3 case is achieved 
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FIG. 7: Unification of the gauge couplings at the two-loop level for central values of low-energy ob- 
servables The SM case with n = 1 is presented by dashed lines. Solid lines correspond to the 

n = 1 scenario with <!>(,, £3 and $ below the GUT scale. Vertical lines mark the relevant scales: Mz, 
M 4i , = M Es = 250 GeV and M GC /t = 3.28 x 10 14 GeV. 

at the two-loop level in the scenario where <£> a is at intermediate and S3 is at the GUT scale. See 
Figs. ^ and El for the n = 1 and n = 2 cases, respectively. As discussed in the text, intermediate 
scale for <£> a would require either small Yukawas for Majorana neutrinos or small C3. In the latter 
case, the novel contributions towards proton decay would be automatically suppressed. The former 
case could be probed if and when the leptoquarks are detected since some of the rare processes 
involving neutrinos would be significantly suppressed compared with the charged lepton ones. We 
note that in the scenario with intermediate $ a the GUT scale grows with the number of light Higgs 
doublets in contrast to the case when E3 is at the intermediate scale. Again, the reason is that 
the B\2 coefficient of <I> a is the same as the appropriate coefficient of the Higgs doublet but its 
contribution to B23 is four times bigger. Thus, the very efficiency of $ a in improving unification 
makes its impact on Mqut rather small. 
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